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ABSTRACT 



Observations of annihilation radiation in the Galaxy 
are briefly reviewed. We summarize astrophysical 
mechanisms leading to positron production and re- 
cent estimates for production rates from nova and su- 
pernova nucleosynthesis in the Galaxy. The physical 
processes involved in the production of annihilation 
radiation in the interstellar medium are described. 
These include positron thermalization, charge ex- 
change, radiative recombination, and direct annihi- 
lation. Calculations of 27 and 37 spectra and the 
positronium (Ps) fraction due to the annihilation 
of positrons in media containing H and He at dif- 
ferent temperatures and ionization states are pre- 
sented. Quenching of Ps by high temperature plas- 
mas or dust could account for differences between 
0.511 MeV and 37 Ps continuum maps. These re- 
sults are presented in the context of the potential of 
INTEGRAL to map sites of annihilation radiation in 
the Galaxy. Positron production by compact objects 
is also considered. 

Key words: positrons; gamma rays; nucleosynthesis; 
supernovae. 



1. INTRODUCTION 



With an omnidirectional flux of « 2.6 xlO -3 ph cm~ 2 
s _1 , the e + -e~ annihilation line at 0.511 MeV is the 
brightest gamma-ray line in the Galaxy. Associated 
with the 27 line is the 37 positronium continuum 
formed by the decay of Positronium (Ps) atoms in 
the 3 Si state. The fraction of positrons that annihi- 
late via the formation of Ps depends on the temper- 
ature, density, composition, and dust content of the 
medium. Consequently observations of 27 and 37 
fluxes and line widths reveal the ISM properties at 
sites where positrons annihilate. Past missions have 
provided limited information on ISM conditions at 
localized sites because they employed non-imaging 
detectors with large fields-of-view (FoVs). Moreover, 
many earlier annihilation line observations used scin- 
tillation detectors rather than germanium spectrom- 
eters, and so were not able to resolve the width of the 



0.511 MeV line. The capability of INTEGRAL to re- 
solve and image the annihilation line — though not 
without considerable technical hurdles — promises 
to revolutionize our view of the Galaxy in the light 
of the e + -e _ annihilation line. 

This review summarizes the status of annihilation 
line astrophysics in the period following the demise of 
the Compton Gamma Ray Observatory and preced- 
ing the launch of INTEGRAL. In Section 2, past ob- 
servations of annihilation radiation in the Galaxy are 
briefly reviewed. Section 3 summarizes the principal 
astrophysical processes that yield positrons, includ- 
ing recent calculations of positron production by no- 
vae, supernovae, cosmic rays, and galactic compact 
objects. Physical processes leading to the 27 line 
and 37 continuum are summarized in Section 4, and 
new calculations of line and continuum spectra and 
37/27 flux ratios in different media are presented in 
Section 5. The review is concluded in Section 6 with 
a discussion of prospects for INTEGRAL to identify 
the origin of positrons that produce annihilation ra- 
diation, including positron production by Galactic 
compact objects. 



SPECTRA OF GALACTIC ANNIHILATION 
RADIATION 



The spectrum of the Galaxy as measured by a soft 7- 
ray telescope is complicated not only by the internal 
and particle backgrounds of the detector, but also 
by variable Galactic 7-ray point sources and the un- 
known distribution of diffuse emissions in the Galaxy. 
These technical difficulties have been dealt with in a 
variety of ways, leading to our present incomplete 
picture of the annihilation emission of the Galaxy. 

The 0.511 MeV line was originally discovered by th e 
Rice University group (Johnson fc Haymes, 1973|), 



though positive identification of t he line required 



the fl ight of a germanium detector ( Lcvcnthal et al _ 
1978| ) which showed that the line, averaged over the 
TFFoV of the Bell-Sandia telescope, has a FWHM 
width AE < 3.2 keV. Other observations of annihila- 
tion radiation have employed broad field-of-view in- 
struments, the largest being the 130° FWHM FoV of 
the Solar Maximum Mission Gamma Ray Spectrom- 
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the galactic CO emission (Share et al., 1990; Harris 
et al., 1990 . T he omnidirectional flux inferred (Pur 



collet al.7l997 ) through analyses of OSSE and earlier 



annihilation line observations is 
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Figure 1. OSSE spectrum of the Galactic Center re- 
gion during Viewing Periods 5 and 16, showing the 
0.511 MeV line and the residual radiation separated 
into the 3"f positronium continuum and background 
e missions. 



ph cm s , consistent with the SMM result. 

Naively assuming that the bulk of the annihilation 
takes place at the distance of the Galactic Center im- 
plies a total Galactic annihilation rate of iV2 7 ~ 4-7r(8 
kpc) 2 (/>27 ~ 2 x 10 43 ph s _1 . Letting / represent 
the fraction of positrons that annihilate via the in- 
termediate formation of Ps, then each e + produced 
2[(1 — /) + (//4)] annihilation photons, because a 
Ps atom decays 1/4 of the time via the 27 state. 
If / = 0.9, as implied by spectral analyses, then 
positron sources in the Milky Way are producing a 
time-averaged rate 



3 x 10 43 e+ s- 1 



The structure in the distribution of annihilation ra- 
diation ca n be determ ined either through m odel- 
dependent (jSkibo. 1993 ) or model- independent ( Pm- 



ccll et al., 19971 ; |Cheng et al, 1997| ; |Milne et al., 199!i| ) 



eter ( Share et al., 1988 ). By showing that the mea- 
sured flu x correlates with the opening a ngle of the 
detector, Lingenfelter fc Ramaty (1989 ) firmly es- 
tablished that there exists a diffuse component of the 
annihilation radiation in the Galaxy. They also ar- 
gued for the existence of a variable, time-dependent 
source of annihilation radiation on the basis of mul- 
tiple observations made with the HEAO-3 and Bell- 
Sandia germanium s pectrometers. Later an alysis of 
the HEAO-3 data ( Mahoney et al., 1994] ) showed 
that the evidence for a variable point source was 
marginal, and no evidence for a variable anni hilation 



mapping techniques. In the former approach, one as 
sumes that the annihilation radiation is described by 
model spatial distributions, for example, those de- 
scribing novae, pulsars, hot plasma, or > 100 MeV 
gamma-ray distributions. One then minimizes \ 2 of 
the observations with respect to various superposi- 
tions of the model distributions. This method is lim- 
ited if the model distributions provide a poor char- 
acterization of features peculiar to the true spatial 
distribution of annihilation radiation. These meth- 
ods were used to argue that a disk component, a 
bulge feature related to an old stellar population, 
and a time- variable point source were required to fit 
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mponont h ao boon o btained through 5. 5 (gw § ^if ble °- 5U MeV Hne observations (fiamatyetal 



e j al,, 1997| ) and 8,5 ([Milne et al,, 1999| ) yeacs-at 



1994) 



observations of the Galactic Center region with the 
OSSE instrument on Compton GRO. 

OSSE has been used to measure the Galactic anni- 
hilation radiation in the greatest detail (Milne et al., 
these proceedings). With a FWHM FoV of 11.4° x 
3.8°, representing about 0.1% of the full sky, OSSE 
can map annihilation radiation through multiple ob- 
servations in different directions. This requires offset 
pointing to subtract background and sophisticated 
mapping techniques to map diffuse emission. An ob- 
servation lasting for two 2-week periods is shown in 
Fig. 1, and illustrates the dominant features in the 
spectrum of a Galactic Center pointing. Shown is the 
spectrum measured during Viewing Periods 5 and 16 
flKinzcr et al., 1999| ). The prominent 0.511 MeV an- 
nihilation line has a flux of ^27 = 1-8 ± 0.19 x 10 -4 
ph cm -2 s . The deduced flux of the 37 continuum 
is </>3 7 = 10.50 ± 0.64,pwr97 x 10~ 4 ph cm- 2 s _1 . 

The omnidirectional flux of 0.511 MeV annihilation 
radiation measured with the Solar Maximum Mis- 
sion is 02 7 = 23tg x 10~ 4 ph cm- 2 s _1 , though the 
precise flux measurement depends on whether the 
underlying distribution of the annihilation radiation 
is assumed to be described by a Galactic Center point 
source or a distribution described by, for example, 



Model-independent mapping techniques, employing 
Singular Value Decomposition or the Maximum En- 
tropy Method algorithms, can be used to reconstruct 
the image without recourse to assumptions about 
und erlying model distri butions. Using these meth- 
ods, Purccll et al. (1997) showed that OSSE observa- 
tions w«re~cIScriEecl~Dy the following emission com- 
ponents: 



• Positive galactic plane: 2.9 ± 0.2, 

• Negative galactic plane: 3.1 ±0.3, 

• Central bulge: 6.0 ± 0.4, 

• Positive latitude enhancement: 2.2 ± 0.2, 

• Mirror region: 0.8 ±0.1 , 

where the fluxes are given in units of 10 ph cm -2 
s _1 . This analysis results in a flux for the full map of 
22.5 ± 0.7 ph cm- 2 s _1 , which is consistent with the 
omnidirectional flux quoted earlier when considera- 
tion of emission in regions not observed with OSSE 
was included. The asymmetry associated with the 
positive latitude enhancement has been interpreted 



Table 1. Positron production processes and mean 
production Lorentz factors 



Processes 


(7) 


N* -> N+c+ 


few 


N + p — ► 7T+ -> C+ 


^60 


N + N' -> N* -» N + c+ 


<L few 


e + B -> 7 + B -> c+c~ 


- 10 6 -10 8 


77 — > e + e~ 


^, few 



( Dormer fc Skibo, 1997 ) in terms of annihilation 
fountain of pair-laden hot plasma driven by a star- 
burst from the Galactic Center region into the halo 
of the Galaxy. 



field near a pulsar polar cap. These positrons can 
escape through magnetic fields that open through 
the light cylinder. The Lorentz factors of positrons 
made in this way are, however, usually so large that 
they will not thermalize on a time scale correspond- 
ing to the age of the Galaxy. The 77 pair production 
process of the fifth row occurs in luminous compact 
regions that may be found in galactic black hole can- 
didates, micro quasars, active galactic nuclei, and in 
gamma-ray bursts (Dermer & Bottcher, these pro- 
ceedings) . Not listed in Table 1 are higher-order pro- 
cesses such as triplet pair production 7+ e — ► e+e + 
+e~, and direct pair production c+ N — > c+ N + 
e + + e _ , which are generally much less important 
because of their small cross sections. 



3.1. Positron Production in Novae 



Later analysis ( Milne ct al., 1999|) has confirmed that 



the positive latitude enhancement has a flux </> 7 = 1- 
2xl0 -4 ph cm -2 s _1 , though there is only marginal 
evidence for a high-latitude enhancement in the 37 
positronium map. This implies that Ps is quenched 
in the fountain, either because the temperature is 
very high or the region is dusty. The first possibility 
can be determined by measuring the line width of 
the annihilation radiation in this direction, and will 
be discussed in more detail in the Section 5. 



There are four principal channels for the production 
of positrons in Novae; these are, 



1. 13 Ni 



.13 



C, r = 598 s, 100% p- 



2 is F ^18 o, r = 109.8 min, 96.9% (3+ 

3. 22 Na ^ 22 Ne, r = 2.60 yr, 90.4% f3+ 

4. 26 A1 ^ 26 Mg, t = 7.2 x 10 5 yr, 82% /?+ and 
7(1.809 MeV), 



POSITRON PROCESSES AND 
PRODUCTION SITES 



As we have seen, sources must produce « 3 x 10 43 
e + s _1 to account for the diffuse annihilation glow 
of the Milky Way. Table 1 lists astrophysical pro- 
cesses that make positrons. The first row represents 
channels where e + are produced through f3 + decay 
of nuclei formed in explosive and hydrostatic nucle- 
osynthesis of novae, supernovae (SNe), Wolf-Rayet 
and Asymptotic Giant Branch stars. The positrons 
produced through this channel have a characteristic 
(3 + decay spectrum with typical kinetic energies of 
~ 1 MeV. The second row represents secondary pro- 
duction interactions between cosmic ray protons and 
ions that produce a characteristic 7r + decay spectrum 
of positrons peaking at ~ 30 MeV. Positrons can also 
be produced in stellar flares through this process. 
The third row represents cosmic ray interactions 
that form ions which subsequently decay through 
the emission of positrons. Depending on whether 
the ions are at rest or are in motion, positrons are 
formed with Lorentz factors ranging from only a few 
to values roughly equal to the Lorentz factors of the 
excited ion. This process can be important in stellar 
flares and in accreting neutron stars and black holes, 
and in regions with enhanced fluxes of low-energy 
cosmic rays. 

The fourth row represents positrons formed by pul- 
sar electromagnetic cascades. Here curvature or syn- 
chrotron photons materialize in the strong magnetic 



where we also list the half-life r and the fraction 
of decays that yield positrons. Nucleosynthesis in 
novae is discussed at this conference by Hernanz; 
here we only note that the first two channels of 
(3 + production are sufficiently prompt that the pro- 
duced positrons will annihilate in t he dense mate- 



rial fo rmed in the nova explosion Hernanz et al 
(1999J ), whereas e + made through the latter two 
channels will be mixed in the expanding envelope and 
can be injected into the ISM. The Galaxy-averaged 
production rate of nov ae from Oxygen/Neon white 
dwarfs is estimated by ( Hernanz ct al., 1999a ) to be 
V 22 = 
V 26 



5 x 10 39 - 2.4 x 10 4U V OA r e/yr e+ s" 1 and 
7 x 10 38 - 4 x 10 40 N O Ne/yr e+ s" 1 for the 
22 Na and 26 A1 channels, respectively. The quantity 
NoNe/yr is the yearly rate of novae throughout the 
Galaxy. |Starrfield et al. (1997| ) give V 22 = 8 x 10 38 

' 26 ^6x 10 39 - 



5 X 10 40 N ON e/yr e+ S 



1 and N'f 
1 for these rates, which are 



2 x 10 41 N ONe/yr e+ s _ 
slightly larger but consistent given the many uncer- 
tainties in the calculation. They also calculate the 
e + production rate from novae on CO white dwarfs, 



6 x 10 37 V 



x l0 37 N co/yr e+ s 



CO/yr e s 

1 . The 



and obtain V 2 2 = 6 x 10 36 

and JVf ^ 7 x 10 36 

yearly rate of novae throughout the Galaxy is esti- 
mated from observations to be « 35 per year. If a 
significant fraction of these novae are produced on 
O/Ne white dwarfs, then novae can provide a non- 
negligible amount of the positrons injected into the 
Galaxy and would contribute to both the disk and 
bulge annihilation radiation. 



3.2. Positron Production from SNe and Stellar 
Winds 



The dominant channels for positron production and 
the associated 7 ray lines formed through hydrostatic 
and explosive nucleosynthesis are 



1. 57 Ni^ 57 Co, r = 36 hr, 40% (3+ and 7(1370 
keV), followed by 57 Co -^ 57 Fe, r = 271 d, 7(122, 
136 keV); 

2. 56 Ni^ 56 Co, t = 6.1 d, 7(150, 750, 812 keV), 
followed by 56 Co ^ 56 Fe, r = 78.8 d, 18% /3+ 
7(847, 1238 keV); 



'3. 44 Ti^ 44 Sc, r S 62 yr, 7 (67.8, 78.4 keV), iol- 
lowed by 44 Sc^ 44 Ca, r = 3.93 yr, 95% /?+ and 
7(1.16 MeV); 

4. 26 A1^ 26 Mg, r = 7.2 x 10 5 yr, 82% /?+ and 
7(1.809 MeV). 

The j3 + decay lifetime from the channel 1 is too short 
to permit the positrons to escape from the ejecta into 
the ISM, so only a prompt 0.511 MeV line could be 
formed, and this would be severely degraded in opti- 
cally thick material. The positron production rates 
into the ISM from Type la and Type II supernovae 
for the remaining channels depend sensitively on the 
percentage r?_2 of positrons that escape int o the ISM. 
For channel 2, Chan fc Lingcnfcltcr (1993 ) find that 



iVl 6 ' 7a = 1.4xl0 4 V2 



MmNi a /c e + s 1 , where M 56 



is the average number of solar masses of synthesized 
56 Fe per SN, and Nj a /c is the mean num ber of Type 
la SNe per ce ntury in the Milky Way. Woosley & 
Weaver (1992| ) calculate that 0.6 < M 56 < 0.9. From 
these rates, it is clear that Type la SNe are capable 
of making an important if not dominant contribu- 
tion to positron product ion, depending cruc ially on 
the escape fraction 77^2 (Milne et al., 1999a). 



Positron production through the 56 Ni— > 56 Co— ► 
channel in Type II SNe is estimated 
Chan fe Lingcnfcltcr (1993|) to be N+' 11 



56 Fe 
by 



xl0 4 V 2 (M 56 /0.08)iV ///c e+ s _1 . Because of the 
thickness of the overlying envelope in Type II SNe, 
the escape fraction is very poorly known and is 
probably <C 1%. At best, this channel could make a 
minor contribution to the total positron production 
rate in the Galaxy. 

The much smaller synthesized mass of 44 Ti com- 
pared with 56 Fe means that the 44 Ti-> 44 Sc-^ 44 Ca 
channel makes a small contribution to the Galac- 
tic positron production rate from Type la SNe, even 
taking into account the longer decay time scale and 
larger positron escape fraction. The production rate 
from Type lb and II SNe depends on the large un- 
certainties in the amount of synthesized mass per 
SNe. In the extreme case, the positron production 
rate from Type la, lb, and II SNe is estimated to be 
j\r 44 ~ 1 n43 



All positrons produced through the 26 A1— > 26 Mg 
chain escape into the ISM due to the long lifetime of 
26 Al. The 1.809 MeV line therefore provides a tracer 
of positron production through channel 4. With 
an omnidi rectional 1.809 M eV flux of « 3 x 10~ 4 
cm~ 2 s~ 4 ( Dichl ct al., 1995J ), this implies an associ- 
ated Galactic positron production rate of ~ 2 x 10 42 
e + s _1 . Type II SNe produce positrons from 26 A1 
at a rate of Nf « 10 41 {M 26 /10~ 5 )N II/C e+ s" 1 , 
where the mass of 26 Al produced per S N ranges from 
0.3 < M 26 /10- 5 < 20 ( prantzos, 1996J ). It is there- 
fore likely that other sources contribute to 26 A1 pro- 
duction, such as massive stellar winds and novae. 
Winds from Wolf-Rayet stars produce positrons at 
the rate of A 26 « 1 x 10 42 



c + s 1 (Meynet et al. 



1997 ). AGB stars (Mowlavi, these proceedings) and 



novae can also make significant contributions. Thus 
the annihilation line map of the Galaxy should dis- 
play a component that traces the 1.809 MeV map, 
although transport of the positrons from the produc- 
tion to the annihilation site could modify the spatial 
correlation, as would be the case for the positive lat- 
i tudc enhancement if d ue to a starburst/superwind 
(permer fc Skibo, 1997[). 



3.3. Other Positron Sources 



Cosmic rays produce positrons through the same 
strong interactions that produce it gamma rays. 
Thus the ir° 7-ray emission centered at sw 70 MeV 
provides a t racer of positron pro duction. COS- 
B analyses (Blocmcn et al., 1984) imply that the 
Galaxy produces ~ 2.5 x lO 4 ^ >100 MeV ph s _1 . 
The 7T° and tt + production cross sections are roughly 
equal, though two 7-ray photons are made from 
each 7T°. Thus the >100 MeV Galactic gamma- 
ray production rate also roughly represents an up- 
per limit to the positron production rate from cos- 
mic ray interactions, after correcting for the emission 
below 100 MeV and noting that bremsstrahlung and 
Compton photons substantially contribute to the dif- 
fuse Galactic gamma-ray background. Large fluxes 
of low-energy cosmic rays would also contribute to 
positron production through (3 + decay of radioactive 
nuclei. COMPTEL uppe r limits to nuclear deexcita- 



tion gamma rays (Blocmcn ct al., 1999) suggest that 
positrons made through this process do not, how- 
ever, significantly contribute to the overall Galactic 
annihilation line flux. 

Although difficult to quantify, compact objects could 
eject large numbers of positrons in to the ISM. Galac 



tic b lack-hole candidate sources (Dcrmcr fc Liang 
1988) ) and microquasars ( Li fc Liang, 1996] ) are the 



two leading candidates for producing localized and 
time-variable sites of annihilation radiation. In- 
dividual sources could eject ~ 10 41 -10 42 e + s _1 
into the ISM. Evidence has been presented from 
Sigma/Granat observations of IE 1740.7-2942 and 
Nova Muscae for episodes of line emission at rs 400 
keV that was attributed to redshifted annihilation 



10 4a e+ s- 1 (Chan fc Lingcnfclter, 1993) 



radiation (For a review, see Harris, 1997), but OSSE 



has not detected variable annihilation radiation or 
annihilation lines from black-hole candidate sources 



such as Cygnus X-l ( Phlips ct al., 1996 ). As already 
noted, electromagnetic cascades in pulsars would 
probably produce positrons with Lorentz factors too 
great to thermalize within the lifetime of the Galaxy. 
Production of positrons by GRBs in the Galaxy 
is discussed in a separate contribution (Dermer & 
Bottcher, these proceedings). 



Table 2. Binary interactions involving positrons. 



Energy Redistribution Positronium Formation 



H 
H 
He 
He 



e+ H* 

e+e- H H 



,+ 



Hc+ 



e e - 
e+ H - 
e+ He 



7 Ps 
H+ Ps 
-> He+ Ps 



Positronium Quenching Annihilation 



4. POSITRON ANNIHILATION PROCESSES 



Early calculation s of positron annihilati on in the ISM 
were treated by |Bussard et al. (1979 ), and in So- 
lar flares by jUranncll ct al. (1976 J^ Updated cross 
sections and processes related to the annihilat ion of 
positrons arc discussed in detail elsewhere ( |Gucs- 
Boum et al., 1991 , 1997 ), so they are reviewed only 



briefly here. Positrons injected into the ISM will 
annihilate directly or thermalize. Thermalization is 
more important than annihilation for positrons in- 
jected at MeV energies, and direct annihilation of 
higher energy positrons produces a line that is invari- 
ably too broadened to be detectable. Positrons ther- 
malize through Coulomb losses with free electrons, 
and through excitation and ionization energy losses 
with bound electrons. When a positron's energy 
reaches a few hundred eV, the energy-loss processes 
compete primarily with the formation of Ps through 
charge exchange in flight, as the rate of direct 
annihilation-in-flight is small. Thermal positrons 
either annihilate directly with free or bound elec- 
trons, or form Ps through radiative recombination 
and charge exchange. Charge exchange with neutral 
hydrogen has a 2430 A threshold, so this process is 
suppressed in sufficiently low-temperature gases. 

Ps is formed in the triplet ( 3 Si) ortho-Ps and singlet 
(^Sq) para-Ps ground states in a 3:1 ratio, with decay 
lifetimes of 1.4 x 10~ 7 s and 1.25 x 10~ 10 s, respec- 
tively. Cascading to the ground state can produce a 
Ly-a line at 6.8 eV and other li nes that are very weak 



but, i n principle, detectab le ( Burdyuzha fc Kauts 
1994 |Wallyn et al., 1996[) . Ps can be quenched if 
the density is ^ 10 13 cm~ 3 , though this is not rel- 
evant for annihilation in the I SM. Annihila t ion of 
thcrmalizcd positr ons on dust (Zurek, 1985; |Gucs- 
30um et al., 1991) is also an important quenching 



process for Ps, leading to a suppression of the 37 
continuum. These calculations are complicated by 
grain properti es, including their sizes , charges, and 
compositions. Guessoum ct al. (1991) find that the 
presence of dust in the warm envelopes around cold 
molecular cloud cores can reduce the Ps fraction, but 
that dust is unlikely to affect the annihilation line 
and continuum properties in the uniform ISM. 



e- 3 Ps - 


-> e e e + 


e + e ^27 


H 3 Ps- 


► H e-e+ 


e+ H -> H+ 2 7 


e~ 3 Ps - 


-+ e~ ^s 


e+ He -> He+ 2 7 


H 3 Ps - 


► H x Ps 


iPs -> 2 7 


H+ Ps - 


-> H e+ 


3 Ps -> 37 


He+ Ps 


-> He e+ 





5. CALCULATIONS OF ANNIHILATION LINE 
AND CONTINUUM SPECTRA 



We have developed a computer code to calculate the 
0.511 MeV line spectrum and relative strength of the 
37 continuum. The code incorporates updated cross 
sections to calculate thermally-averaged rates for all 
relevant binary interactions involving positrons, and 
uses a Monte Carlo technique to simulate the ther- 
malization and Ps production of ^MeV positrons. 
We include charge exchange with neutral and singly- 
ionized helium, which have generally been ignored in 
astrophysical applications. Processes involving dust 
and H2 are not yet included. 

For the formation of Ps, we include radiative re- 
combination with free electrons and charge exchange 
with hydrogen and neutral He. For Ps quenching, we 
include breakup of Ps by neutral hydrogen, ortho-Ps 
quenching via free electrons and via electrons bound 
in H, inverse charge exchange with H and He, and 
ionization of ortho-Ps atom by free electrons. This 
last reaction turns out only to be important in ion- 
ized gases at densities £ 10 13 cm -3 , and is relevant 
to positron annihilation at the Sun. 

For the Monte Carlo determination of the fraction of 
positrons forming Ps before thermalization and the 
resulting Ps energy distribution, we include He and 
we use an updated experimentally-determined cross 
section for charge exchange in hydrogen. We find 
that this new cross section results in higher Ps pro- 
duction than previously calculated, and the resulting 
annihilation line width is broader. 

The Ps formation, quenching and energy loss pro- 
cesses along with the positron annihilation processes 
treated in our calculations are shown in Table 2. The 
important diagnostic spectral features of annihila- 
tion radiation are the ratio Q?,~jlQi~i and the FWHM 
width AE of the line. These two quantities depend 



Table 3. Temperatures and densities, in units of 
cmT 3 , of different phases of the ISM used in calcu- 
lations of annihilation radiation 





Cold 


Warm 


Warm 


Hot 




Neutral 


Neutral 


Ionized 




T(K) 


80 


8000 


8000 


4.5 x 10 5 


nrn 


38 


0.31 


0.06 


0.0 


n e 


0.0 


0.055 


0.17 


3.5 x 10~ 3 


hhi 


0.0 


0.055 


0.17 


2.9 x 10~ 3 


IlHel 


3.8 


0.036 


0.023 


0.0 


nHoII 


0.0 


0.0 


0.0 


0.0 


nHoIII 


0.0 


0.0 


0.0 


0.6 x 10~ 3 



upon the conditions of the medium where the annihi- 
lation occurs. The medium can be characterized by 
the temperature T, the total (neutral + ionized) hy- 
drogen density njj, and the ionization fractions X H + , 
X He + and X He ++. We use the non-thermal and 
thermally-averaged rates discussed above to solve a 
system of continuity equations for pair equilibrium 
and provide rates for positron annihilation from all 
of the various processes. These rates then directly 
give Q3-y/Q2-y The total line spectrum is the sum of 
several components, some having Gaussian shapes, 
others not. Hence, the total line shape will not, in 
general, be Gaussian. 

Fig. 2 shows calculations of the annihilation line and 
Ps continuum for phases of the ISM with tempera- 
tures and densities given by Table 3. Fig. 3 shows 
calculations of the same quantities for high temper- 
ature, fully ionized media. The quantities Q3-y/Q2-y, 
giving the relative number of photons in the 37 Ps 
continuum to the 0.511 MeV line, and the FWHM 
line width AE, from Figs. 2 and 3, are plotted in 
Fig. 4. In the limit of a fully neutral medium, all 
annihilation takes place through the formation of Ps. 
Consequently Qd, 1 /Q2-y — ► 9/2, because 3/4 of the 
annihilations take place via ortho-Ps, giving 3 pho- 
tons per annihilation, whereas 1/4 of the annihila- 
tions take place via para-Ps, giving 2 photons per 
annihilation. Annihilation occurs before thermaliza- 
tion in a fully neutral medium, so the calculated 9 
keV line width for the neutral medium is broader 
than for partially ionized media where the positrons 
generally thermalize before annihilating. 

As the temperature increases, the relative num- 
ber of 37 continuum to 27 line photons decreases. 
This is because the Ps formation processes involv- 
ing charge exchange with neutrals become less im- 
portant at higher temperatures due to the fewer 
available neutral or partially ionized ions. The line 
width increases at high temperatures due to thermal 
broadening, and approaches the limit AEfwhm 
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Figure 2. Calculations of annihilation line and Ps 
continuum from the injection of Me V positrons into 
regions with different temperatures and ionization 
fractions. 
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1)97), where Tq is the temperature in units of 10° K. 
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Fig. 4 quantitatively shows how the 37 Ps contin- 
uum is quenched in high temperature plasmas. Con- 
sequently Ps formation can be quenched in high tem- 
perature gas. This property can account for the 
marginal detection or absence of a high-latitude en- 
hancement in the Ps maps of the Galactic Center 
region as compared with the observation of a high- 
l atitudc enhancem ent in the 0.511 MeV line maps 
(Milne ct al., 199E). If this interpretation is correct, 
then a broadened annihilation line would be detected 
in the direction of the high-latitude enhancement 
with a line width AEfwhm ~ 10 keV and with an 
intensity ~ 1-2 x 10~ 4 ph c m' 2 s -1 . For broad FoV 
instruments such as TGRS ( Harris ct al., 199S ) , this 



broadened line is difficult to resolve from the narrow, 
bright annihilation line emitted by the disk and bulge 
components. The imaging capability of INTEGRAL 
will be able to determine whether the interpretation 
of the high latitude enhancement as a fountain of 
hot gas venting into the halo of the Galaxy is correct 
by mapping the intensity and line width in this di- 
rection, provided that the positive latitude enhance- 
ment is sufficiently localized. Further work will be 
required to determine if Ps can be quenched by dust 
in the Galactic Center region, which could also lead 
to suppression of the 37 Ps continuum. 



Figure 3. Calculations of annihilation line and Ps 
continuum from the injection of Me V positrons into 
high temperature plasma. 
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Figure 4- Calculations of 37/27 ratios and FWHM 
line widths in different phases of the ISM. 



6. SUMMARY AND CONCLUSIONS 



We have summarized our knowledge of astronomical 
sources of positrons and the processes that lead to 
the formation of annihilation radiation in the ISM. 
Because of the large uncertainties in nucleosynthesis 
calculations, the relative contributions of novae, SNe, 
and Wolf-Rayet winds remains unclear. Positrons 
from cosmic rays, compact objects, and GRBs could 
also make a significant contribution to positron pro- 
duction in the Galaxy. An important tool to reveal 
the dominant sources of positrons will be to use 7- 
ray tracers to establish associated positron injection 
rates. COMPTEL and INTEGRAL maps of the 26 Al 
1.809 line and EGRET maps of tt° gamma-ray emis- 
sion should trace some fraction of the annihilation 
glow of the Galaxy, with the residual attributed to 
other processes. 

One of the lingering uncertainties in positron as- 
trophysics is the question of whether positrons are 
formed by accreting compact objects in the Galaxy. 
SPI on INTEGRAL will be well suited to examine 
this question, not only by searching for annihilation 
line transients, but also by looking for annihilation 
hot spots that are coincident with compact object 
emission in the Galaxy. SPFs 511 keV sensitivity 
for on- axis point sources dur ing a 10 6 s observatio n 
is 2-3 x 10~ 5 ph cm" 2 s" 1 flVcdrcnne et al., 1999J ). 
Spectral analysis may show that annihilation pho- 
tons from directions coincident with known compact 
objects have a different line width or 37/27 ratio 
than the surrounding emission, suggesting a point 
source origin. For sufficiently bright sources, SPI 
could provide sub-degree imaging that could be com- 



pared with IBIS maps of point sources. Searches for 
variations of the Ps fraction and line width will re- 
veal if there are unusual sites of positron production, 
and whether the emission from the positive latitude 
enhancement arises from annihilation in a hot super- 
wind or at compact objects. 
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